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The Fermi LAT collaboration has recently presented constraints on the gamma-ray signal from 
annihilating dark matter using separate analyses of a number of dwarf spheroidal galaxies. Since 
the expected annihilation signal has the same physical properties regardless of the target (except 
for a normalization scale), it is possible to enhance the constraining power using a combined 
analysis, the initial results of which will be presented here. 
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1. Introduction 

The Fermi Gamma-ray Space Telescope was launched on June 11, 2008. Its main instrument, 
the Large Area Telescope (Fermi-LAT), consists of 16 identical modules in a 4x4 array with each 
module comprised of a tracker for direction determination and a calorimeter for energy measure- 
ments [1]. The Fermi-LAT observes the entire sky every ~3 hours (2 orbits) with a field of view 
covering ~2.4 sr and a sensitive energy range extending from 20MeV to >300GeV. These features 
makes the Fermi-LAT a great instrument for dark matter (DM) searches. 

The leading DM candidate is a weakly interacting massive particle (WIMP), and this is the 
DM candidate we are primarily focusing on with the Fermi-LAT. The gamma-ray flux from self- 
annihilating WIMPs can be expressed as ^imp{E, y) = J{v) ^ '^'^''{E) [2], where <^^{E) is the 
"particle physics factor" described by 

^PP 1 <av> ydN^ 

and J{y) is the "astrophysical factor", or J-factor, described by 

JiY)= f dliY)pHK¥))- (1-2) 

Jl.o.s. 

Here (av) is the WIMP annihilation cross section times relative velocity, mwimp is the WIMP mass, 
p (r) is the dark matter density distribution, and ^ j- f is the gamma ray spectrum generated per 
WIMP annihilation where the sum is over final states / with branching ratio Bf. The particle 
physics factor has two spectral features: the continuum feature and the line feature (which is often 
referred to as the "smoking gun"). 

Dwarf spheroidal galaxies (dSphs) are DM dominated systems since they have a high mass to 
light ratio. Many dSphs are closer than 100 kpc to the galactic center, and they have low background 
since most dSphs are expected to be free from other astrophysical gamma-ray sources and they have 
a small gas content. This makes them a good target for gamma-ray DM searches with the Fermi- 
LAT. However, the photon statistics for single dSphs are expected to be very low. The Fermi-LAT 
collaboration has recently presented results from a DM search in a number of dSphs [3]. Since 
the DM spectra are the same for all astrophysical sources, a combined analysis would improve the 
statistics. 

In this proceeding we present the preliminary results of a combined likelihood analysis of eight 
dSphs. 

2. Analysis 

The combined likelihood used in this work is described by 

N 

L{av,mwiMp;obs) = Y[Li{av,mwiMP,C,bi;obs), (2.1) 

i 

where av and mwimp are the common DM signal parameters (velocity averaged cross-section 
and mass), C ai^e constants (e.g. branching fraction in our case), and b are individual parameters 
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Name 



1 



b 



19 GeV^ 



cm 



Bootes I 358.08 

Coma Berenices 241.9 

Draco 86.37 

Fornax 237.1 

Sculptor 287.15 

Sextans 243.4 

Ursa Mayor II 152.46 

Ursa Minor 104.95 



69.62 

83.6 

34.72 

-65.7 

-83.16 

42.2 

37.44 

44.80 




Table 1: dSphs analyzed in this work. All positions and J-factors of the dSphs are taken from [3]. 

(isotropic and galactic diffuse background normalization, normalization of nearby sources). The 
main advantages of the combined likelihood are that the analysis can be individually optimized and 
that combined limits are more robust under individual background fluctuations and under individual 
astrophysical modelling uncertainties than individual limits. 

We have observed the eight dSphs listed in Table 1, which is the same subset analyzed in [3]. 
In this preliminary analysis we have used 21 month data from 2008-08-04 to 2010-05-12. We have 
used the diffuse event class which only contains the events with the highest gamma-like confidence, 
and we have chosen events ranging from 200MeV to lOOGeV. We used the Fermi-LAT instrument 
responce function P6_V3_DIFFUSE. Our region of interest (ROI) is a region of 10 degrees radius 
centered on dSph location. Standard cuts removing Earth albedo photons have been made. The 
dSphs are modeled as DM point sources using the DMFit package [4] where we only consider 
100% annihilation into the bb annihilation channel so far since it represents the spectral shape 
of a large fraction of the possible final states. The background is modeled according to Fermi- 
LAT recommendations [5], and nearby soures are modeled according to the first year point source 
catalog [6]. We perform a binned analysis to use both energy and spatial information. The data 
selection and analysis were performed using the Fermi-LAT analysis package, ScienceTools [7], 
and the upper limits are obtained using profile likelihood as implemented in the MINUIT processor 
MINOS [8]. 

In this proceeding work we have not yet taken systematic uncertainties into account. The 
uncertainties mainly arise from the J-factors (about 20 to 200% as can be seen in Table 1), and the 
effective area of the LAT (5 to 20% and energy dependent [9]). However, some initial tests show 
that the combined results are less sensitive to J-factor changes than individual limits. 

3. Preliminary results 

In Fig. 1 we present constraints on the WIMP annihilation cross-section for eight dSphs, 
both individual and combined limits, evaluated at mwjMp = 50, 100, 150, 300, 600 and 1000 GeV. 
No systematic uncertainties have been included and the results ai^e preliminary. The combined 
limits approach the parameter space below the thermal WIMP cross-section for low WIMP masses. 
Limits down to 5 GeV were presented at the conference, but due to an extrapolation problem for 
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Figure 1: Upper limits on WIMP annihilation cross section for annihilation into 100% bb, av evaluated at 
m-wiMP = 50, 100, 150, 300, 600 and 1000 GeV. The expected thermal WIMP cross-section is plotted as a 
reference. The limits are 84% two-sided. 



masses below 10 GeV only limits for mass points above 10 GeV are reliable. An update for low 
masses is forthcoming. 

4. Discussion and outlook 

A stacking analysis of 8 dSphs has been presented using a combined likelihood approach. 
Limits improve with respect to the the most stringent of the 8 individual limits, depending on 
WIMP mass. Some tests on consistency under choice of ROI, fit range and binning have been 
performed. The results presented here are preliminary. The inclusion of systematic uncertainties is 
in progress. More dSphs will be added if they are sufficiently nearby, if relatively accurate estimates 
of their DM distribution can be obtained, and if they are situated at sufficiently high galactic latitude 
to avoid galactic foreground. The analysis will be updated for more recent data, more annihilation 
channels will be studied, and a paper is in preparation within the Fermi-LAT collaboration. 
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